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Abstract. Using the most recent updated physics, cal- 
ibrated solar models have been computed with the new 
thermonuclear reaction rates of NACRE, the recently 
available European compilation. Comparisons with mod- 
els computed with the reaction rates of Caughlan & Fowler 
( |1988p and of Adelberger et al. ( p.998[ ) are made for global 
structure, expected neutrinos fluxes, chemical composition 
and sound speed profiles, hclioscismological properties of 
p-modes and g-modes. 

Key words: Physical data and processes: nuclear reac- 
tions, nucleosynthesis, abundances - Sun: evolution - Sun: 
interior - Sun: Oscillations 



1. Introduction 

Precise solar models have been constructed over the past 
three and a half decades ( see, e.g 

1998b| ) 



Bahcall et al. 1963 



> . Bahcall 1993, Bahcall et al 



The refinement have 
accelerated in the past decade (see, e.g. Bahcall and Ul- 
ri 



ich 1984 Bahcall & Pinsonneault 1992b). Since ten years 
many stellar and solar models have been computed using 
the thermonuclear reaction rates of Caughlan & Fowler 
( |1988| , hereafter C88) a popular compilation but not op- 
timized for solar conditions. More recently some authors 
(e.g. Bahcall et al. 1995; Reiter et al. 1995| ; Chaboyer et 
al. 1995|; Berthomieu et al. 1995; Christensen-Dalsgaard et 
al. [1996 ) employed the improved thermon uclear reaction 
rates adopted by Bahcall & Pinsonneault ( 1992a ) for the 
calculation of accurate solar models. Meanwhile several 
groups of nuclear physicists have undergone other compi- 
lations of updated thermonuclear reaction rates of astro- 
physic al int erest. A year ago the compilation of Adelberger 
et al. ( 1998 , hereafter A98) has been published. The origi- 
nal motivation of this compilation is to assess the state of 
the nuclear physics important to the solar neutrino prob- 
lem. The incidences on solar models of these new reaction 
rates have been a nalyzed by se veral groups (see e.g . Bah - 
call et al. |l998b| ; Brun et al. |l998| ; Morel et al. |1998| ). 
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More recently the European Nuclear Astrophysics Com- 
pilation of REaction rates (Angulo et al. 1999| , NACRE, 
hereafter N99) has been completed and opened to free 
access. The driving motivation of this last work, coordi- 
nated by the Institut d'Astronomie et d'Astrophysique of 
the Universite Libre de Bruxelles, is the build-up of well 
documented and evaluated sets of experimental data or 
theoretical predictions of astrophysical interest. To have a 
direct idea of the degree of reliability of any reaction rate, 
the authors have published, either a very convenient plot 
of the available measurements of cross section S-factors 
with respect to energy, or a table showing the range of 
the various parameters needed for cross section evaluation, 
e.g. the resonance parameters. Moreover, the accuracy of 
each analytical fit is indicated. This new compilation gives 
besides the adopted reaction rate 1Z, its lower and upper 
limits 1Z\ and 1Z U . A solar neutrino analysis based on pre- 
liminary NACRE d ata for the PP reactions has bee n done 
by Castellani et al. (|1997| ). Recently Arnould et al. ( [1999D 
have used the N99 reaction rates to compute abundance 
predictions in non hydrogen and helium burning. They 
convincingly show that large spreads in the abundances 
predictions for several nuclides may result not only from 
a change in temperature, but also from nuclear physics 
uncertainties. 

We are now in the fortunate position of having two 
precise and independent determinations of the best nu- 
clear fusion data, namely A98 and N99. In order to illus- 
trate the effects, on the standard solar model, of nuclear 
fusion rates on various astronomical quantities, including 
neutrino fluxes and helioseismology frequencies, we com- 
pare the model results calculated with the best current 
data from A98 and N99, with results obtained using early 
estimates of fusion rates of C88. Those differences are not 
very large. Nevertheless they modify the energy balance, 
the stratification, the chemical composition and the neu- 
trino generation in the core. 

Let us first recall the main constraints known nowa- 
days on solar models. The helioseismological constraints 
relevant to the core are the small p-mode frequency dif- 
ferences Sv 02 and Svi 3 and the not yet observed spec- 
trum of gravity modes. Other signatures of changes of 



2 



P. Morel, B. Pichon, J. Provost & G. Berthomieu: Solar Models and NACRE thermonuclear reaction rates 



Table 1. Respectively for the C88, A98 and N99 compilations, the S-factors at zero energy S(0) (MeV barn) and 
S'(0) (barn) are given for each nuclear reaction of solar interest, but 7 Be(e~, iA- Bc 7) 7 Li. The last column gives A1Z, 
the global uncertainty of N99's rates computed, for T 6 = 15, according to Eq. [TJ. The S(0) factor of 12 C(p, 7) 13 N is 
estimated from the N99 plot. 
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thermonuclear reaction rates will be the sound velocity 
profile which is known from inversion of helioseismic data 
between R <; 0.1 i?© and R ^ 0.9 Rq, and also the ra- 
dius of the base of the solar convection zone which is pre- 
cisely located. The amount of observed photosphcric de- 
pletions of lithium and beryllium which are often ascribed 
to transport phenomena beneath the convection zone are 
also sensitive to changes of the thermonuclear reaction 
rates in their low energy regime. An another constraint 
also connected to nuclear reaction rates is the isotopic ra- 
tio 3 He / 4 He measured at present day at the solar surface 
which is sensitive to the pre-main sequence deuteron burn- 
ing and to initial isotopic ratios 2 H / X H and 3 He / 4 He 
of cosmological interest. 

The paper is organized as follows: in Sec. || in the low 
energy range, for the nuclear reactions of interest for solar 
modeling, we summarize the main differences between N99 
and A98 with respect to C88. The physics used in the 
models is described Sec. E| In Sec. |] we report results of 
comparisons between calibrated solar models computed 
with N99, A98 and C88, finally we conclude in Sec. |. 



2. Comparison of thermonuclear reaction rates 
from N99, A98 and C88 compilations 

The most important nuclear reactions relevant of solar 
modeling are, for the PP chains (Clayton 1968| ; Bahcall 



1989 Table 3.1 and 3.3): 



^(^/J+i/pp) 2 !!, 2 H(p, 7 ) 3 He, 3 He( 3 He, 2p) 4 He, 
3 He(a, 7 ) 7 Be, 7 Be(e", V7 Bcl ) 7 U, 7 Li(p,a) 4 He, 
7 Be(p, 7 ) 8 B*(/3+!/s B ) 8 Be(a) 4 He, 

and for the CNO bi-cycle: 

12 C(p, 7) 13 N(/?+^3 N ) 13 C, 13 C(p, 7 ) 14 N, 

14 N(>, 7) 15 0(/3+M5 ) 15 N, 15 N(p, 7) 16 0, 15 N(p, a) 12 C, 

16 0(p, 7 ) 17 F(/3+M7 F ) 17 0, 17 0(p,a) 14 N. 

Owing to their low termination and small contribution 
to energetic and nucleosynthesis, despite their interest 
for neutrino generation, we do not explicitly take into 
account in the nuclear network 1 H(pe~, t/ pcp ) 2 H and 
3 He(p, e + fhep) 4 He the so-called pep and hep reactions. 
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Fig. 1. Left panel: relative differences between the reaction rates between N99 and C88 for the relevant solar temper- 
ature range 0.5 < T 6 < 19. for the PP reactions: ^(p, /3 + z/ pp ) 2 H (heavy, full), 2 H(p,7) 3 He (dashed), 3 He( 3 He, 2p) 4 He 
(dot-dash-dot-dash), 3 He(a, 7) 7 Be (dotted), 7 Li(p, a) 4 He (dash-dot-dot-dot) and 7 Be(p, 7) 8 B* (thin, full). Right panel: 
same between A98 and C88. 



Nevertheless we compute the number of f pep neutrino gen- 
erated using the equation (3.17) of the Bahcall's (1989) 
reference text book. 

The changes between the reaction rates of N99, A98 
and C88 are extensively commented in Adelberger et al. 
( |1998| ) and Angulo et al. ( p,999[ ). As a matter of illustra- 



3 He( 3 He, 2p) 4 He: For the most energetic reaction of 



tions, for the three compilations and for each PP and CNO 
reaction - but the electronic capture 7 Be(e~, iA- Bc 7) 7 Li - 
Table gives the S-factors at zero energy and the under- 
lying global uncertainty on the rate AH: 



ATI 



1. 



(1) 



Tlx 

estimated for Tq — 15; T§ is the temperature in MK, TZ\ 
(resp. 7l u ) stands for lower (resp. upper) limit of N99 up- 
dated reactions. For our thermonuclear reaction network 
the contributions of resonances to the astrophysical reac- 
tion rates are negligible in the solar range of temperatures, 
therefore the values of S(0) and, if any, S'(0) presented 
here are pertinent. For sake of briefly we omit to repro- 
duce the known S" (0) values. Figure [l] (resp. Fig. ||) com- 
pares the relative differences between the adopted rates 
of N99 (resp. A98) and C88 for the temperature range 
0.5 <T 6 < 19. 

We next briefly recall the main changes in the rates of 
A98 and N99, with respect to those of C88 which is the 
oldest and, up to nowadays, the most used and complete. 

2 H(p, 7) 3 He: Among all reactions of PP chains and 
CNO bi-cycle, it is the rate of this PPI reaction which is 
the most badly known. The reaction rate is so fast that it is 
only involved by the pre-main sequence deuteron burning. 
Owing to the lower value adopted for the S-factors at zero 
energy, the rate of the reaction which synthesizes 3 He is 
about —24% lower in N99 than in C88. This reaction is 
not updated in A98, for the calculations with A98 we used 
the value adopted in C88. 



the PP chains N99 (resp. A98) adopts values smaller by 
about -6% (resp. -2%) than C88 for the S-factors. As 
consequences of the calibration process, for the models us- 
ing either N99 or A98, more 1 H nuclear fuel will be burnt 
in order to reach, at present day, the observed luminos- 
ity and effective temperature. Therefore these models will 
have cores with larger temperature, helium content, den- 
sity and sound velocity than models computed with C88; 
then, at first sight, their predicted total neutrino fluxes 
are expected to be larger. This effect will be enhanced for 
the models computed with N99 since the rates of the two 
reactions 1 H(p, /3+f pp ) 2 H and 2 H(p, 7) 3 He are smaller in 
N99 than in C88. 



Figure p. 26 of Angulo et al. ( 1999| ) gives the impres- 
sion that, at low energy, the Junker at al.'s ( 1998[ ) recent 
measurements are avoided by the interpolation formulae 
adopted (see also Fig. 2 of Adelberger et al. 1998j ). These 
recent data should lead to an increase of S(0) so, to an 
enhancement of the efficiency of the reaction and then, 
owing to the calibration process, to a decrease of solar 
neutrino fluxes. 

7 Li(p, q) 4 Hc : The S(0) value adopted by N99 differs 
from C88 between +15% and +7%. In the core its curate 
rate is irrelevant due, first to its strong rate (~ 10~ 5 year, 



e.g. Bahcall 1989 Table 3.2) and, second to the tiny mass 
fraction of 7 Li ~ 210~ 15 . Beneath the convection zone 
the burning of 7 Li will be more efficient with N99 than 
with C88, leading to an increase of the lithium depletion 
at the solar surface at present day. As suggested in C88, 
the 7 Li burning is slightly enhanced by few percents by 
the neighbor reaction 7 Li(p, 7) 8 Be(a) 4 He which has been 
added to our nuclear network. This reaction is not updated 
in A98, for the calculations with A98 we shall use the value 
adopted in C88. 
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Fig. 2. The same as Fig. | for the CNO reactions: 12 C(p, 7 ) 13 N (heavy, full), 13 C(p, 7 ) 14 N (dashed), 14 N(p, 7 ) 15 
(dot-dash-dot-dash), 15 N(p, 7 ) 16 (dotted), 15 N(p,a) 12 C (dash-dot-dot-dot), 16 0(p, 7 ) 17 F (thin, full), 17 0(p,a) 14 N 
(heavy, dashed). 



7 Be(e , ^7 Bc7 ) 7 Li: N99 deals only with charged par- 
ticle induced reactions involving nuclei, and therefore the 
7 Be electron capture rate is not updated. In the calcu- 
lations with N99 we shall used the value given by A98. 
Beneath Tq = 1 only an upper limit is given in C88 for 
the 7 Be electron capture. The adopted rate of A98 differs 
from the rate of C88 by more than +50% at low temper- 
ature; for Tq ~ 15, i.e. in the solar core, the rates of C88 
and A98 are of same order. 

7 Be(p, 7 ) 8 B*: This reaction controls the efficiency of 
the important source of ^sb, the so-called boron solar neu- 
trino. The adopted values for the S-factors at zero energy 
are slightly larger in N99 than in A98, but still smaller 
than in C88. With respect to C88, everything else equal, 
one can expect that the neutrino flux from boron will be 
reduced for the solar models computed with A98 and N99. 

13 C(p, 7 ) 14 N: The values of S-factors at zero energy 
adopted by N99 and A98 are magnified by a factor of 
about two with respect to their previous values in C88; 
as a result the rates are increased by +30% and +15% 
respectively. These large differences will not have any no- 
ticeable incidence on the global structure of the core since 
the energy generated by the CNO bi-cycle is only ^ 2% 
of the total nuclear energy. 

14 N(p, 7 ) 15 0: The rate of the most important reaction 
for the computation of energy generation and neutrino 
fluxes created by the CNO bi-cycle is known with a large 
uncertainty. The three compilations adopt about the same 
values for the S-factors at zero energy. Figure shows 
small differences between the rates. This is due to different 
interpolation formulas which slightly differ since there is 
no measurement at low energy (see the convincing figure 
p. 58 of Angulo et al. p99| )). 

15 N(p, 7 ) 16 0: For the reaction which governs the effi- 
ciency of the NO-part of the CNO bi-cycle, N99, A98 and 
C88 adopt the use of S-factors at zero energy obtained by 



Rolfs & Rodney ( 1974 ). Due to differences in the interpo- 
lation formula Fig. |2j reveals enhanced rates of +15% in 
N99 with respect to C88 or A98. 

16 0(p, 7 ) 17 F: At low energy the reaction which con- 
trols the generation of w? F , the so-called fluorine solar 
neutrino, is based on data with large experimental errors. 
The adopted rate has the largest uncertainty among the 
CNO reactions. Though Table [l] gives for the three com- 
pilations about the same values for the S-factors, Fig. || 
shows large differences for the rates resulting from dif- 
ferent analytical formulatio ns. B eyond Tq ~ 10, N99 and 
C88 are close (Angulo et al. |l999|) . The difference of -50% 
between A98 and C88 results of the used of the stan- 
dard formulation of the non-resonant reaction rate with 



S-factors (Fowler et al. 1967 ). 

17 0(p, a) 14 N: N99 and A98 use different a nalyt ical fits 
based on the measurements of Landre et al. (1989). They 



differ by —30%. With the discovery of a resonance at low 
energy (Landre et al. loc. cit.) the analytical fit of C88 be- 
came in error by more than two order of magnitude. For 
the models computed with C88 we have used the rates 
derived from the Landre's et al. analytical fit, as recom- 
mended by A98. 

Summary. With respect to C88 many reaction rates, prin- 
cipally 3 He( 3 He, 2p) 4 He, are lowered in N99 and also, but 
in a less extend, in A98. One can expect that this will lead 
to calibrated solar models with central cores with larger 
temperature. For the reactions of PP chains, with respect 
to C88, other important changes connected to the observ- 
able neutrino fluxes are the rates of the electronic cap- 
ture 7 Be(e~, JA- Bo7 ) 7 Li which is significantly diminished 
for Tq ^ 7 in A98 and, for N99, the decrease of the rate 
of 7 Be(p, 7 ) 8 B*. With respect to C88, the changes in N99 
and A98 of the reaction rates of the CNO bi-cycle are not 
large enough to modify significantly the solar model. 
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Fig. 3. Relative differences in opacity k (full), heavy element Z (dashed), density p (dot-dash-dot-dash), temperature 
T (dotted) and helium content Y (dash-dot-dot-dash) for the calibrated models A98 (left) and C88 (right), with 
respect to model N99. 
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Fig. 4. Relative differences of abundances with respect to model N99, for models A98 (left) and C88 (right) for PP 
species: 1 U x 100 (full), 3 He x 10 (dashed), 4 He x 100 (dash-dot-dash), 7 Li (dotted), 7 Be (dash-dot-dot-dash). 



3. The solar models 

Basically the physics of the models is the same as in Morel 



et al. (1997) 



Calibration of models. Each evolution is initialized with a 
homogeneous zero-age pre-main-sequence model in quasi- 
static gravitational contraction with the temperature at 
center T c ~ 0.5 MK, i.e. close to the onset of the deuteron 
burning. The models are calibrated within a relative ac- 
curacy better than 10 -4 by adjusting: the ratio l/H p of 
the mixing- length to the pressure scale height, the ini- 
tial mass fraction A; of hydrogen and the initial mass 
fraction {ZjX\ of heavy element to hydrogen in order 
that, at present day, the solar models ha ve the luminosity 
L Q = 3.846 1 33 erg s^ 1 (Guenther et al. |l992|) , the radius 



Rq — 6.9599 10 10 cm (Guenther ct al. loc.cit.) and the 
mass fraction of heavy el ement to hydrogen (Z/X)q — 
0.0245 (Grevesse & Noels 1995 ). We used a time of evo- 
lution < cv = 4600 My, an intermediate value between the 
mete oritic age i 0m = 4530 ± 40 My of the SunQ (Guen- 
ther 1989 ) and its helioseismic value t Q h = 4660 ± 100 My 
derived by Dziembowski et al. ( 1998 ). The zero age main- 
sequence (ZAMS) is defined as the time where nuclear re- 
actions dominate gravitation as the primary energy source 
by more than 50% (Guenther et al. loc. tit.). The mass of 
the Sun is assumed to be M Q = 1.9891 10 33 g (Cohen & 



Taylor 1986) 



1 Here tQ m is referenced with respect to ZAMS which occurs 
36 ± 10 My (Guenter loc. cit.) after th e form ation of meteorites 
4566 ± 5 My from now (Bahcall et al. |l995|). 
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Nuclear and diffusion network. The general nuclear net- 
work we used contains the following species : 1 H, 2 H, 3 He, 
4 He, 7 Li, 7 Be, 9 Be, 12 C, 13 C, 14 N, 15 N, 16 0, 17 and 
Ex; Ex is an "Extra" fictitious mean non-CNO heavy el- 
ement with atomic mass 28 and charge 13 (Ex ~ 28 Al) 

which complements the mixture i.e., Xgx = 1 — zJ =ih ^ 
with Xi as the mass fraction of the species labeled with 



H, 



17 



O. With respect to time, due to microscopic 



diffusion processes, the abundances of heavy elements are 
enhanced toward the center; Ex mimics that enhancement 
for the non CNO metals which contribute to changes of 
Z , then to opacity variations but neither to nuclear en- 
ergy generation nor to nucleosynthesis. To compute the 
depletion of 9 Be, we have added, to the nuclear network 
given Sec. ||, the most efficient reactions of 9 Be burning: 
9 Be(p, d)2 4 He and 9 Be(a,n) 12 C. Th e life time of the neu- 
tron, namely 888s (Barnctt et al. 1996 ), is smaller by 
more than thirteen orders of magnitude than the evolu- 
tionary time scale of the Sun's main-sequence pp reaction. 
Therefore, for the calculations, the last reaction is rewrit- 
ten 9 Be(a, e~pp9Bc) 12 C. The weak screening of Salpeter 
( |1954| ) is used, it is a very good approximate of the exact 
solution of the Schrodinger equation for the fundamental 
pp reaction (Bahcall et al. 1998a). 



The protosolar initial isotopic ratios (in number) for 
hydrogen and helium are respectively taken as 2 H/ 1 H = 
3.01 10~ 5 , 3 He/ 4 He= 1.1 IO -4 (Gautier & Morel ^99l\) . 
The initial ratios between the heavy elements within Z are 
set to their photosphcric present day values, namely (in 
number) C: 0.24551, N: 0.06458 and O: 0.51295 (Grevesse 



& Noels 1993) then, for the complement Ex: 0.17696. The 



initial isotopic ratios are deri ved f rom the abundances of 
nuclides (Anders & Grevesse [l989l) : 13 C/ 12 C = 1.11 10" 3 , 
15 N/ 14 N= 4.25 IO" 3 , 17 0/ 16 0= 3.8110- 4 . We have used 



the meteoritic values (Grevesse & Sauval 1998j ) for the 
initial abundances in dex, (H = 12), of Li and Be: 



Li 



= 3.31 ±0.04, 



Be 



H 



= 1.42 ±0.04. 



For the calculations of depletions, the lithium is assumed 
to be in its most abundant isotope 7 Li form, so it is with 
beryllium assumed to be 9 Be. Neither the meteoritic abun- 
dance nor the nuclide isotopic ratio of 7 Be are known, due 
to numerical constraints the protosolar abundance of 7 Be 
was somehow arbitrarily taken to a very low, but non zero 
value, namely [ 7 Be/ 1 H] = —3.58 dex. The initial abun- 
dance of each isotope is derived from isotopic ratios and 
initial values of X = X H + 2 H, Y = 3 He ± 4 He and Z/X 
as inferred by the calibration process in order to fulfill the 
basic relationship X + Y + Z = 1. 

Microscopic diffusion is described by the simplified for- 
malism of Michaud & Proffitt (1993) with each of the 
heavy elements as a trace element. 



Equation of state, opacities, convection and atmosphere. 
We have used the OPAL equation of state (Rogers et al. 
1996) and opacities (Iglesias & Rogers 1996) for the solar 
mixture of Grevesse & Noels ( |1993[ ) complemented, at low 
temperatures, respectively by the MHD equation of state 
(Dappen 1996) and Alexander & Ferguson ( [19941 ) opaci- 
ties. The interpolations of opacities are made with the v9 
birational spline package of G. Houdek (Houdek & Rogl 



1996; Houdek 1998) 



In the convection zones the temperature gradient is 
computed according to the standard mixing-length theory. 
The mixing-length is defined as I = aH p , where H p is the 
pressure scale height. The convection zones are mixed via 
a strong turbulent diffusion coefficient, which produces a 
homogeneous composition. 
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The atmosphere is restored using a T(r) law derived 
from an atmosphere model of the Sun computed by van't 
Veer ( |1998[ ) with the Kurucz's fll99ip ATLAS 12 package. 
The connection with the envelope is made at the Rosse- 
land optical depth Tb = 20 (Morel et al. 1994), where 



the diffusion approximation for radiative transfer becomes 
valid. A smooth connection of the gradients is insured be- 
tween the uppermost layers of the envelope and the opti- 
cally thick convective part of the atmosphere. The radius 
R+ of any model is taken at the optical depth ~ 0.54 
where T(t±) = T c s; the mass of the star M*, is defined as 
the mass enclosed in the sphere of radius i?* . The exter- 
nal boundary is located at the optical depth r cxt = 10~ 4 , 
where the density is fixed to its value in the atmosphere 
model p(T ext ) = 3.55 10~ 9 gcm~ 3 , that corresponds about 
to the temperature minimum in the solar chromosphere. 



Numerics. The mode ls ha ve been computed using the 
CESAM code (Morel 1997 ). The numerical schemes are 
fully implicit and their accuracy is first order for the time 
and third order for the space. For numerical performance 
and algorithmic constraints the analytical expressions of 
reaction rates are tabulated with respect to temperature 



for the range 0.5 < Tq < 20 and interpolated with a 
relative accuracy better than 10 -5 . Each evolution needs 
about 90 models. Typically 600 mass shell are used along 
the evolution, it increases up to 2100 for the models used 
in seismological analysis. 



p-mode and g-mode oscillation calculations. The frequen- 
cies of linear, adiabatic, global acoustic modes of the so- 
lar models have been computed for degrees £ = to 
£ = 150 and have been compared to the observations. 
The characteristic low degree p-mode frequency differ- 
ences Av nt f — v n j — Vn-i.t+2 for £ — and £ = 1, which 
provide information on the properties of the solar core, 
have been fitted by linear regressions with respect to n: 

&v n ,i = 6u n ,i + Si(n - no), n = 21, £ = 0, 1, 

both for the observations and the theoretical frequencies. 
For the gravity modes which have not yet been observed, 
we give the characteristic asymptotic spacing period Pq 
according to Provost & Berthomieu ( |1986| ). 

4. Comparison of models 

Table || gives the global properties of models and Fig. || 
exhibits the profiles, with respect to radius, of the most 
important variables for the internal structure namely, den- 
sity, temperature, opacity, helium and heavy element con- 
tents. 

4-1. Chemical composition 

The changes in chemical composition directly result from 
changes of thermonuclear reaction rates but also, in a more 



intricate way, from changes in microscopic diffusion coeffi- 
cients which are sensitive to the temperature and density, 
and to chemical composition, pressure, temperature and 
density gradients. 



Changes at the surface and in the envelope. For the three 
models N99, A98 and C88, Table | shows that the ex- 
pected photospheric abundances of helium are slightly re- 
duced and remain compatible with the range of observed 
values. As known (Basu |l997b ), the amount of photo- 
spheric "observed" helium derived from inversion of he- 
lioscismic data is more sensitive to the equation of state 
than the amount of photospheric "predicted" helium de- 
rived from calibrated solar models. Indeed we have cali- 
brated a solar modelQ using C88 thermonuclear reaction 
rates and the MHD (Dappen 1996) equation of state in- 
stead of OPAL and obtained a photospheric helium con- 
tent Y s = 0.246 which is the value derived from inversion 



using the MHD equation of state (Basu & Antia 1995). 

Though the 7 Li surface depletion is increased with the 
use of the enhanced rate of 7 Li(p, a) 4 He adopted by N99, 
the predicted abundance is still very far from the ob- 
served value. The differences are likely consequences of 
the lack, in standard solar models, of mixing generated by 
the shear at the level of the tachocline (see e.g. Gough et 
al. 1996, Brun et al. 1999). It smoothes the chemical com- 



position gradients and reduces the microscopic diffusion 
efficiency (Basu 1997a) immediately beneath the convec- 
tion zone. According to the new observations (Grevesse 
& Sauval 1998), the predicted photospheric depletion of 
beryllium is tiny. The predictions for the surface isotopic 
ratios ( 3 He/ 4 He) s by the three models are all within the 
interval of accuracy given by the observations. 



Changes in the core. The solar core is the innermost part 
where the nuclear energy generation is efficient. It extends 
from the center to about R c ~ 0.4i? Q slightly beyond 
the 3 He peak located around 0.3i? Q . Owing to the less 
efficiency of PP reactions, see Sec. 0, the temperature, the 
density, the amount of helium and the sound velocity at 
center of calibrated models N99 and A98, are larger than 
in C88. As expected, Fig. || shows almost symmetrical 
profiles for the differences for 1 H and 4 He. Owing to the 
larger efficiency of the pp reaction in C88, larger values are 
obtained for the relative difference C88 minus N99 than 
for A98 minus N99. The typical features for the relative 
differences of abundances of 3 He are consequences of the 
smaller reaction rates of N99 with respect to A98 and C88 
of 3 He( 3 He, 2p) 4 He. Beneath the 3 He peak, owing to the 
increase of the temperature, the amount of 3 He is smaller 
in C88 and A98 than in N99; it is the reverse beyond 
the peak. Though the same rate prevails in N99 and A98 
for the 7 Be electronic capture, there is a non zero value 



Not analyzed here for sake of briefness. 
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Table 2. Comparison of global characteristics of solar models computed with the thermonuclear reaction rates of 
respectively Angulo et al. ( |1999| , N99), Adelberger et al. ( |199S| , A98) and Caughlan & Fowler ( |1988[ , C88) compilations. 
The first four rows give the initial values of Yi (resp. Z\) the mass fraction of helium (resp. heavy element), {Z/X)i the 
mass fraction of heavy element to hydrogen, a the mixing-length parameter. 7 Li sz {resp. 9 Be sz ) and 7 Li s [resp. 9 Be s ) 
respectively are the surface depletions in dex (H = 12) of 7 Li {resp. 9 Be) at zero-age main sequence and at present 
day. All other item are for present day. The surface isotopic ratio ( 3 He/ 4 He) s is in unit of 10~ 4 ; Y s , Z s respectively 
are the surface mass fraction of helium and heavy element; Rcz is the radius, in solar units, at the bottom of the 
convection zone; T c , p c , Y c and Z c are the central values respectively of, the temperature in units of 10 7 K, the density 
in gem -3 , the mass fraction, of helium and heavy element. Svo2 and 8v\^ are the values, in /iHz, of the frequency 
differences between the radial p-modes of degree I — — 2 and I = 1 — 3. Po is the characteristic spacing period of 
g-modes in minutes. 



N99 A98 C88 Observed values 





0.2723 


0.2726 


0.2729 


Zi 


0.0197 


0.0197 


0.0196 


{z/xy 


0.0278 


0.0278 


0.0277 


a 


1.924 


1.931 


1.941 




2.26 


2.38 


2.37 


9 Be sz 


1.42 


1.42 


1.42 


7 Li s 


2.18 


2.30 


2.29 


9 Be s 


1.35 


1.35 


1.353 


( 3 He/ 4 He) a 


4.34 


4.32 


4.32 


Y s 


0.2436 


0.2442 


0.2447 


z s 


0.0181 


0.0181 


0.0181 


Rcz 


0.7138 


0.7132 


0.7124 


T c 


1.573 


1.570 


1.566 


Pc 


153.8 


153.0 


151.9 


Y c 


0.6418 


0.6420 


0.6409 


Zc 


0.0210 


0.0210 


0.0210 


5VQ2 


9.21 


9.18 


9.16 




16.10 


16.06 


16.03 


Po 


35.13 


35.23 


35.42 



1.10 ±0.10, Grevesse & 
1.40 ±0.09, Grevesse & 
4.4 ± 0.4, Bodmer et al 
0.232 



0.249, Basu (1997b) 



Sauval ( 1998 ) 
S auva l ( J1998I ) 
(|l995|) 



0.713 ± 0.001, Basu & Antia (1995) 



9.002 ± 0.044 - 9.014 ± 0.042, from LOI/GOLF (see text) 
15.884 ± 0.034 - 15.711 ± 0.071, from LOI/GOLF (see text) 



for relative difference between the 7 Be profiles of models 
N99 and A98, resulting of differences between the rates 
adopted for 3 He(a,7) 7 Be. 

Figure || exhibits large differences in the abundances of 
16 O for A98 and C88 with respect to N99 despite the fact 
that the rates of the reactions of 16 burning are close. 
In fact, these differences result from changes of rates of 
15 N burning which creates 16 0. For 12 C and 14 N, around 
0.18-R©, effects of changes of nuclear reaction rates are 
magnified by the large gradients of that species. There, 
Fig. |3| reveals, on Z profiles, small bumps due to the mag- 
nification by large gradients of variations in chemical com- 
position caused by the changes of thermonuclear reaction 
rates. 



4-2. Thickness of the convection zone. 



For radius R il 0.4i? Q , i.e. in the envelope, Fig. |^ shows 
that the opacity profiles are close within ±0.4% for mod- 
els N99 and A98. The thickness of the convection zone is 
about the same in N99 and A98, and close, within the er- 
ror bars (see Table |^), to the observed value. It is slightly 
larger for model C88. That difference is due to the increase 
of the radiative temperature gradient resulting from the 
higher value of the opacity. The differences of temperature 
between N99 and A98 being small, the changes in opacity 
are mainly due to the variations of density. The relative 
opacity differences amount to ±1% between C88 and N99. 



P. Morel, B. Pichon, J. Provost & G. Berthomieu: Solar Models and NACRE thermonuclear reaction rates 



9 



d 



02 

d 



/ 



0.05 



0.1 

R/R 



0.15 



0.2 



0.05 



0.1 

R/R G 



0.15 



0.2 



Fig. 6. Left panel: relative differences of neutrino rates generation, between models A98 
full), i/7 Bc (dashed), va B (dot-dash-dot-dash), ws N (dotted), ws (dash-dot-dot-dash). Right 
models C88 and N99 



for v pp (thin, 



and N99, 
panel: the same between 



Table 3. The same as Table || for the neutrino fluxes. Us- 
ing the comprehensive symbols, defined Sec. ^, z^ pp , " pcp , 
^ 7 Bej ^Bi t/l3 Ni Vl5 Oi Vl7 F are, at earth level, the number 
cm~ 2 s _1 of neutrinos of each kind. 4>Q a and <&ci m SNU 
and <&Ka m events day -1 , are the expected fluxes for the 
three neutrino experiments namely gallium, chlorine and 
Kamiokande (see text). The observed values for $Gaj ^Ci 
and <£>Ka are respectively from Hampel et al. (1999), Davis 
( |1994| ) and Fukuda et al. ( [19961 ). 



N99 



A98 



C88 



Observed values 



"pp 


5.91 10 10 


5.92 10 10 


5.94 10 10 


Z^pcp 


1.40 10 8 


1.40 10 s 


1.48 10 8 


"'Be 


4.90 10 9 


4.80 10 9 


4.71 10 9 


U8 B 


5.68 10 6 


5.17 10 6 


6.19 10 6 


"13 N 


5.73 10 8 


5.7710 s 


5.34 10 s 


"ISO 


4.96 10 8 


5.01 10 s 


4.5710 s 


Un F 


6.41 10 6 


3.15 10 6 


5.74 10 6 


$Ga 


130.1 


128.4 


129.8 


<I>C1 


8.31 


7.71 


8.82 


$Ka 


0.61 


0.55 


0.66 



77.75 ± 6.2 
2.55 ± 0.25 
0.29 ± 0.02 



+4.3 
4.7 



The hep flux which may be important (e.g. Bahcall & 
Krastev 1998c, Fiorentini et al. 199S) in the neutrino spec- 



trum measurements by the SuperKamiokande, SNO and 
Icarus experiments is not listed. With respect to C88, due 
to hotter core, za- Bo and CNO neutrino fluxes are enhanced 
in A98 and N99 and, as expected, z/ pp is slightly reduced. 

Despite larger temperatures in the core we obtained, 
for the models A98 and N99 with respect to the model 
C88, the expected decreases of the vs B boron neu- 
trino fluxes owing to their reduced rate of the reaction 
7 Bc(p, 7) 8 B*. The introduction of N99 reaction rates rel- 
atively to A98 induces an increase of +10% of ve. B . The 
effect is significant on the flux measured by the chlorine 
and Kamiokande experiments. Note that Table ^ reveals 
that the neutrino fluxes at earth level for A98 are very 
similar than those given Table 1 in Bahcall et al. (1998b). 
They differ only by few percent for vn-p owing to the large 
abundance of 16 resulting from the great efficiency of 15 N 
burning which about double the fraction of termination of 
the NO part of the CNO bi-cycle. Despite the fact that the 
two stellar evolution programs are entirely independent of 
each other, when the same nuclear reaction rates are used 
(A98), the most important neutrino fluxes ^ pp , " pop , "?Bc, 
^8 B , M3N! and "i5Q all agree to better than 1%. 



4-3. Neutrinos 

Table || gives the predicted neutrino fluxes at earth level 
and the expected fluxes for the three neutrino experiments 



namely, chlorine (e.g. Davis 1994), gallium (e.g. Hampel et 
al. 1999) and Kamiokande (e.g. Fukuda et al. 1996[ ), com- 
puted according to Berthomieu et al. (1993). The gallium 



and chlorine absorption c ross s ections have been ta ken re - 
spectively from Bahcall ( 1997 ) and Bahcall et al. ( 1996| ). 



4-4- Seismological comparison 

The seismic properties of the solar model are mainly re- 
lated to the profile of sound-speed (resp. Brunt- Valssala, 
frequency) as far as p-modes (resp. g-modes) are con- 
cerned. Figure [?] shows that the models N99 and, to 
a lesser extend A98, compared to C88, have a larger 
sound speed in the central core below 0.3 solar radius 
by +0.2% (resp. +0.1%), and a smaller one by -0.1% 
(resp. —0.05%) just below the convection zone. Table || 



10 



P. Morel, B. Pichon, J. Provost & G. Berthomieu: Solar Models and NACRE thermonuclear reaction rates 



1 




o 

X 1 












CD 

I 2 - 

i x 


\ : - ...... ^ J\\ 


A 


O N 
\ 

■a 




1 

Q 
U 

1 

o 

X 
EM 

1 



Fig. 7. Relative frequency differences of sound speed 
(full) and square of Brunt- Vaissala frequency divided by 
10 (dashed) between N99 (thick) and C88, and A98 (thin) 
and C88. 
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Fig. 9. Relative difference in sound velocity between the 
Sun and the calibrated models N99 (full), A98 (dashed) 
and C88 (dot-dash-dot-dash). 
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Fig. 8. Frequency differences (in /iHz) in the low fre- 
quency range between N99 and C88 for modes of degree 
£ = 0, 1, 2. Open symbols denote g-modes and, full sym- 
bols, f-modes and p-modes. 



shows a small increase of small low degree differences 61*02 
and <5^i3, as defined Sec. || in relation with the difference 
of sound speed in the solar core. There the relative dif- 
ferences in the Brunt- Vaissala frequency between models 
N99 and C88 are larger by a few percents, i.e. one order 
of magnitude larger than the sound speed differences. The 
increase of the Brunt- Vaissala frequency leads to a —1% 
smaller value of P , the characteristic spacing period of 
g-modes. The differences between A98 and C88 are three 
times smaller. This change in Brunt- Vaissala frequency 
influences much the low frequency modes for frequency 
less than 1 mHz, i.e. the low radial order p-modes, the 
f- and g-modes. Consequently the frequency differences of 
the low degree p-modes between models N99 and C88 vary 
from —0.1 /iHz to —0.25 /iHz when the frequency increases 
from 1 mHz to 5 mHz, with a minimum value of —0.5 /iHz 



around 2 mHz. The normalized frequencies differences for 
p-modes of degree £ — 3 to I = 150 are negative and 
change by less than —1 /iHz in the observed range. As ex- 
pected, the change of nuclear reaction rates do not modify 
the frequency of oscillation of degree larger than 70. 

Figure |8| shows the frequencies differences in low fre- 
quency range between N99 and A98 and Table [| gives 
the frequencies of g- and p-modes for £ — to I = 2 
in the same frequency range. In the low frequency range 
400 //Hz - lmHz, it appears that the p-mode frequencies 
are changed by less than 0.5 /iHz between N99 and C88, 
i.e. by less than +0.1%, with an effect larger at lower fre- 
quency. Below 200 /iHz, the oscillations are gravity modes 
with an asymptotic behavior, and the relative period dif- 
ferences are almost proportional to Pq given in Table | (see 



Provost et al. 



199J for details). Between 200 to 400 /iHz, 
or f- 



the oscillations are gravity modes, or f- and pl-modes 
and they are more influenced by the change of the Brunt- 
Valssala frequency in the solar core, induced by changes 
in nuclear reactions, except the pl-modes for £ — and 
£ = 1. Frequencies shifts are much larger, of the order of 
1 to 1.5 /iHz, when the frequency varies from 200 /iHz to 
400 /tHz, i.e. about 1%, when comparing the models N99 
and C88. 

In the range 0.1 Rq Ss R 0.9 Rq where the inver- 
sions of the helioseismic data are reliable, the sound speed 
of the three models has been compared with the seismic 
sound speed experimental results of Turck-Chieze et al. 

shows that the relative differences are 
. The discrepancy between the Sun and 



( |1997D . Figure | 
below a few 10" 
the models is larger for model N99 with sound speed too 
small just below the convection zone and too large in the 
core. Table || shows that it is the same for the quantities 
5i>02 and 5i>i3 of the models compared to the correspond- 
ing observed values 8v n ^ derived from GOLF (Grec et al. 
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Table 4. Frequency of g-modes and p-modes in the range 
100 fiHz-2 mHz. T is the type of the mode labelled by the 
radial order. 



1 


C88 


N99 


T 


I C88 


N99 


T 





257.82 


258.31 


pl 


2 101.88 


102.64 


glO 




404.32 


404.79 


P 2 


111.07 


111.88 


g9 




535.98 


536.30 


P 3 


121.91 


122.80 


g8 




680.67 


680.71 


P 4 


134.87 


135.82 


g7 




825.56 


825.52 


P 5 


150.51 


151.55 


g6 




972.95 


972.82 


P 6 


169.69 


170.83 


g5 




1118.36 


1118.15 


p7 


193.24 


194.49 


g4 




1263.78 


1263.60 


P 8 


221.23 


222.57 


g3 




1407.86 


1407.58 


P 9 


255.22 


256.67 


g2 




1548.78 


1548.51 


plO 


295.44 


296.93 


gl 




1687.12 


1686.81 


pll 


354.60 


356.06 


f 




1822.51 


1822.14 


pl2 


383.36 


384.60 


pl 




1957.79 


1957.44 


pl3 


514.35 


514.58 


p2 










664.33 


664.33 


P 3 


1 


108.55 


109.31 


g5 


811.77 


811.75 


P 4 




127.04 


127.93 


g4 


959.86 


959.75 


P 5 




152.45 


153.50 


g3 


1105.18 


1105.01 


P 6 




190.56 


191.85 


8'2 


1250.78 


1250.60 


p7 




261.62 


263.07 


gl 


1394.73 


1394.46 


P 8 




284.77 


285.27 


pl 


1536.06 


1535.80 


P 9 




448.34 


448.48 


p2 


1674.78 


1674.44 


plO 




596.90 


596.95 


p3 


1810.37 


1810.01 


pll 




746.65 


746.66 


P 4 


1946.00 


1945.63 


pl2 




893.71 


893.60 


P 5 










1039.63 


1039.52 


P 6 










1185.68 


1185.47 


p7 










1329.80 


1329.57 


P 8 










1473.13 


1472.88 


P 9 










1612.86 


1612.53 


plO 










1749.57 


1749.25 


pll 










1885.34 


1884.94 


pl2 









1997]) and VIRGO/LOI (Frohlich et al. |1997| ) observations 
on SoHO. 



5. Discussion and conclusions 

We have compared the structure, the neutrino fluxes, the 
chemical composition profiles and the helioseismological 
properties of calibrated standard solar models computed 
with the adopted nuclear reaction rates of the European 



compilation NACRE (Angulo et al. pL999[ ) with those of 
calibrated solar models computed with the nuclear reac- 
tion rates of Caughlan & Fowler ( 1988] ) and Adelberger et 
al. ( |1998| ). 

Roughly speaking, the thermonuclear reaction rates of 
PP chains adopted by NACRE and, but in less extend, 
by Adelberger et al., are slightly less efficient than those 
adopted by Caughlan & Fowler. The calibration generates 
models with cores of larger temperature, density, helium 



content and sound speed with the concomitant increase 
of the neutrino fluxes, except for v pp and z/s B ; for this 
last one, the decrease is due to the smaller rate of the 
reaction 7 Be(p, 7) 8 B*. Thus the predicted neutrino fluxes 
are reduced for the chlorine and Kamiokande experiments, 
but almost unchanged for gallium. For Kamiokande and 
chlorine, N99 predicts intermediate values between A98 
and C88. 

The introduction of the NACRE thermonuclear rates 
increases the discrepancy between predicted and observed 
sound velocity profiles between the Sun and the models, 
both below the convection zone and in the solar core. 
These relative differences, though at the level of a few 
thousandths, are smaller for the model computed with 
the reaction rates of Caughlan & Fowler, the increase is 
- +0.5% for C98 and - +1% for N99. The radius at the 
base of the solar convection zone is in good agreement 
with the observed value for all models. 

Though NACRE adopts an enhanced rate for the reac- 
tion of lithium burning 7 Li(p, a) 4 He, the predicted deple- 
tion of photospheric lithium remains too small to fit the 
observed value. 

The differences between calibrated solar models com- 
puted with the adopted thermonuclear reaction rates of 
the two new compilations are rather small. It is not really 
possible to make any choice between them. From an in- 
crease of the accuracy of the observed p-mode frequencies 
and from a hopefully detection of g-modes one can expect 
to improve our knowledge on the stratification of the so- 
lar core in the goal to validate, in the low energy regime, 
the thermonuclear reaction rates and their concomitant 
neutrino generation. 

Thanks are due to the NACRE 's work which also pro- 
vides these estimates of uncertainties on the adopted rates. 
These new features, now available to the user, are impor- 
tant to constraint the solar model. We are investigating 
this point in a work in progress. 
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